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0>-» ' Flow equations for an 0(Af)-syninietric effective potential are discussed and solved 

for the finite temperature case. The model is investigated at the critical point and 
critical exponents for various A'^ are calculated. 

(N. 

,_^ i 1 Introduction 

>: 

\^ , The A^-component scalar field theory with an 0(iV)-symnietry often serves as 

'nI" ■ a prototype for symmetry restoration investigations at finite temperature. The 

'^ , manifestation of universality of critical phenomena enables the applicability of 

^*"^ ' the 0(A^)-symmetric model to a wide class of very different physical systems 

Qs , around the critical temperature. 

0^ I Wilczek and Pisarski have argued that the two flavor chiral transition is 

(— I . either of second order or of first order depending on the C/a(1) anomaly. In case 

pi i' that the transition is of second order Wilczek and Rajagopal have shown that 

JL i two flavor QCD belongs to the same iiniversality class as the three dimensional 

fl_j ' four component Heisenberg magnet El. Thus for A^ = 4 this theory is used as 

r^ . an effective model for the chiral phase transition in two flavor QCD enabling 

^ ' deeper insights into the chiral phase transition in QCD provided that the chiral 

order parameter remains close to thermal equilibrium through the transition. 
On the lattice it is still not clear whether QCD for two flavors and the 0(4)- 
model lie in the same universality class □. 

In a second order transition, the system is at an infrared flxed point of the 
renormalization group. This means that the physics is scale invariant and the 
flow equations show a scaling solution. Due to a diverging correlation length 
in the chiral limit the order parameter fluctuates on all length scales and the 
correlation functions exhibit a power law behavior with critical exponents. 

The article is organized as follows: In Sect. we summarize the concept 
of the renormalization group (RG) method. In Sect. 0we present some results 
for the 0(4)-model at finite temperature. Finally Sect. H is dedicated to the 
critical behavior of the system at exactly the critical temperature. 



- 1 — I 

X 



2 Concept of the Renormalization Group Method 

RG flow equations predict the behavior of a given theory in different momen- 
tum regimes. The main ingredient of the RG approach is the integration of 
irrelevant short-distance modes in order to derive a low-energy effective the- 
ory in the infrared in analogy to a discrete block-spin transformation on the 
lattice. For this purpose one introduces an infrared scale k that separates the 
fast-fluctuating short-distance modes from the slowly- varying modes. This idea 
- pioneered by Wilson and Kadanoff in the early seventies B - can be system- 
atically incorporated e.g. in a momentum regularization in a transparent way. 
It results in an effective action parameterized by the averaged blocked field at 
the scale k. Thus the full one-particle irreducible (IPI) Feynman graphs or 
vertex functions are generated by the renormalized effective action in the limit 
/c — > and the fc-dependent effective action provides a smooth interpolation 
between the bare action defined at the UV scale A, where no fluctuations are 
considered, and the renormalized effective action in the IR. The flow pattern 
of the given theory is obtained by studying an inflnitesimal change of the IR 
scale k in the effective fc-dependent action and is governed by differential flow 
equations. 

Of course, the integration step cannot be performed in an exact way and 
one has to resort to some approximation such as the loop expansion. However 
a sharp momentum cutoff is in conflict with important underlying symmetries 
of the considered theory. The task is then to implement both the UV and the 
IR cutoff scales in a symmetry-conserving way, which can be accomplished by 
an operator cutoff regularization in Schwinger's proper-time representation. 

The one-loop contribution to the effective action in general yields a non- 
local UV diverging logarithm, which can be rewritten in a proper-time repre- 
sentation and leads to 

oo 
1/A2 

where the primes on the potential V denote differentiation w.r.t. </>. Here, the 
UV (IR) divergences appear for t ^ (oo). We modify the above expres- 
sion by introducing a regulating a priori unknown smearing function fk (r) in 
the proper-time integrand. Differentiating the resulting fc-dependent effective 
action w.r.t. the scale fc yields the flow equations. It can be shown that all 
universal results in_the infrared are independent of the choice of this regulating 
smearing function Q 



3 Finite Temperature 

It is straightforward to generalize the above approach to finite temperature. 
Within the imaginary time (Matsubara) approach the four-dimensional the- 
ory at T = is conveniently matched to the effectively three-dimensional 
behavior at the critical temperature by a replacement of the four-dimensional 
momentum integration in Eq. (nh with a three-dimensional integration and a 
Matsubara summation over the corresponding bosonic frequencies. Thus at fi- 
nite temperature the scale k serves as a generalized IR cutoff for a combination 
of the three-dimensional momenta and Matsubara frequencies. 

To solve the flow equations numerically we start the evolution in the broken 
phase deep in the UV region fc = A and use a tree-level parameterization of 
the potential with two initial values. These initial values are fixed at T = 
and are constant in a relatively large temperature region. 
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Figure 1: Left panel: The potential for various temperatures. Right panel: The temperature- 
dependent order parameter in the chiral limit. 

In the fc -^ limit the temperature dependence of the potential V{(1>) 
(left panel of Fig. ^ and the order parameter (j)k^o{T) (right panel) signal a 
second-order phase transition. Our calculations coincide with chiral perturba- 
tion theory up to T ~ 45 MeV. Around Tc we obtain a scaling behavior of the 
order parameter with a critical exponent /3 oc 0.4 (cf. right panel). 



4 Critical Behavior 

At the critical temperature only properly rescaled quantities asymptotically 
exhibit scaling and the evolution of the system is purely three-dimensional d. 
This is the well-known dimensional reduction phenomenon. The rescaled di- 
mensionless minimum K{t) of the potential approaches a constant (the fixed 



point K*) as the dimensionless flow "time" t = ln(/s/A) tends to — cx) in the 
infrared (see left panel of Fig. 0). For a starting value near the critical value 
Kcr at the ultraviolet scale i = the evolution moves towards either the spon- 
taneously broken k 7^ or the symmetric k = phase in the infrared. Due to 
the rescaling, which in the infrared tends to a constant value or zero during 
the fc-evolution, the dimensionless minimum k — po/fc^ diverges for /c — > in 
the broken phase. The time the system spends on this scaling solution can 
be rendered arbitrarily long by appropriate fine tuning of the initial values at 
t = 0. 
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Figure 2: Left panel: The minimum k, r) and coupling A as function of t. Right panel; The 
phase diagram for N = 1 (see text for details). 



Besides of the trivial high-temperature and low-temperature fixed points 
the 0(7V)-symmetric model exhibits a nontrivial mixed fixed point inherent in 
the flow equations. The region around this fixed point is displayed in right 
panel of Fig. where the rescaled dimensionless quantities n and A are plotted 
for iV = 1 and for different initial values during the evolution toward the 
infrared. 

In general there are two relevant physical parameters for the 0(A^)-model 
which must be adjusted to bring the system to the critical fixed point accord- 
ing to universality class arguments. Since we here work in the chiral limit (no 
external sources or masses) only one relevant eigenvalue from the linearized 
renormalization group equations is left (the temperature). This can also be 
seen in Fig. 0. The quantity k is the relevant variable because repeated renor- 
malization group iterations (which are the discrete analog to the continuous 
evolution with respect to i — > —00) drive the variable away from the fixed 
point value while A is the irrelevant variable being iterated toward the fixed 
point if the initial values are chosen sufficiently close to it. Thus there is a one- 



dimensional curve of points attracted to the fixed point which is the so-caUed 
critical surface for this model. 

One nicely sees that this critical line separates both phases. When we 
choose initial values near the critical line the system spends a long time near the 
critical point. E.g. starting with k > Kcr and A finite the evolution ultimately 
tends to the zero-temperature fixed point corresponding to very large values of 
K. On the other hand starting below the critical line the evolution of k tends 
towards the high-temperature fixed point corresponding to small k's. Only for 
exactly k = Kcr at i = the renormalization group trajectory flows into the 
critical mixed fixed point, which means that the long distance behavior at the 
critical point is the same as that of the fixed point. 

Here the predictive power of this nonperturbative approach becomes vis- 
ible: During the evolution near the scaling solution the system loses memory 
of the initial starting value in the UV and the effective three-dimensional dy- 
namics near the transition is completely determined by the fixed point being 
independent of the details of the considered microscopic interaction at short 
distances. Exactly at a second-order phase transition the potential should be 
described by a fc-independent solution. The equation of state drives the po- 
tential away from the critical temperature. As a result, after the potential has 
evolved away from the scaling solution, its shape is independent of the choice 
of the initial values and exhibits universal behavior. 

In the vicinity of the critical point we obtain a scaling behavior of the 
system governed by critical exponents which parameterize the singular behav- 
ior of the free energy around the phase transition. Out of the six critical 
exponents for the 0(iV)-model only two are independent due to four scaling 
relations which we have verified by explicit calculations of a, /3, 7 and v in 
three dimensions for various N (see tab. p. 

The critical exponent /3 parameterizes the behavior of the spontaneous 
magnetization or the order parameter for the broken phase (ka > Kcr) in the 
vicinity of Tc- The difference ka — Kcr is a measure of the distance from the 
phase transition irrespective of any given value of Aa. If ka is interpreted as a 
function of temperature this difference is oc (T ~Tc), i.e. Kcr defines the critical 
temperature in three dimensions. 

All calculated critical exponents nicely converge to the large-N values. In 
addition for iV = 1 we have found 77 = 0.0439 and 6 = 4.748. 

The structure of the regularized RG-improved flow equations depends on 
the used a priori unknown cutoff functions fk- Employing the same class of the 
cutoff function as in ref.u we have calculated critical exponents for several cut- 
off functions with more and more monomials on a grid for the full potential. A 
small systematic decrease in the values for the critical exponents with the num- 
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0.3475 


-0.085 


1.39 


N=3 
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0.79 
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0.4555 
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0.4965 
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1.986 
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1.0 


0.5 


-1.0 


2.0 



Table 1: Critical exponents for different N compared with the large-N result. 

ber of included monomials of the smearing function is observed. Restricting to 
less monomials in the cutoff function /fc seems to accelerate the IR convergence 
at the critical point. This is visible in the evolution of the minimum of the 
potential which stays longer on the critical trajectory in the vicinity of the 
scaling solution. For more monomials the evolution of the minimum escapes 
faster from the critical trajectory to the symmetric as well as the broken phase. 

To conclude we have presented a powerful nonperturbative method based 
on the RG approach applied to the 0(A'^)-model at finite temperature, which, 
in fact, is consistent with pertinent lattice simulations. It remains to be seen in 
how far the critical 0(4)-behavior is realized in full two flavor QCD simulations. 
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